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ABSTRACT 



The work of Lewis and Walecka on the structure of the giant reso- 
nance in Carbon 12 and Oxygen 16 has seemed to indicate that the 
single -particle, shell-model theory gives better results than collective 
models --in particular the Goldhaber - Teller and Steinwedel - Jensen 
models. The transverse form factor as a function of momentum 
transfer is observed experimentally^ fall to a minimum and rise 
again in the range from zero to 0. 6 f . This dip is predicted by the 
Brown particle -hole theory model but is not predicted by either of these 
collective models. The purpose of this paper is to develop a collective 
model which could possibly give this dip in the form factor. Further 
work is necessary to determine if this model can actually do so. A 
semiclas sical treatment is used to find the matrix elements for the 
assumed model. The transverse form factor is then calculated for two 
nuclear charge density functions -- constant density and the Fermi 



distribution. 
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1. Introduction. 



Electron scattering is a very useful tool for the study of nuclear 
structure, since the cross-section for electron scattering from a 
nucleus can be related to the basic geometric details of the structure 
of the nucleus. The incoming electron is assumed to interact with the 
nucleus by means of the electromagnetic field. This interaction is 
known, whereas the nuclear interaction (strong force) is not. 

The parameters of electron scattering are the initial momentum k 
and initial energy E of the incoming electron, which are controlled by 
the experimenter; and the final momentum k 7 and the final energy E 7 , 
which are measured. The cross-section can be found as a function of 
the momentum transfer q - k - k from the electron to the nucleus 
and of the excitation energy, which can be found from E, E 7 , the 
scattering angle, and considerations of the nuclear recoil. One can then 
vary q for a fixed excitation energy (i. e. , for an excitation of the 
nucleus to a particular level) , and get the cross-section as a function 
of q (where q = | q | ) for that level. This cross-section is related to 
the Fourier transform of the charge and current densities of the nucleus. 
Thus, we can get information about the charge and current densities of 
the nucleus from electron scattering. 

A so-called "dipole state 11 is a feature of the electron scattering 
cross-section for most nuclei. The dipole state is a state strongly 
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excited by photons of the appropriate energy (usually 10 - Z5 Mev), It 
shows up as the M giant dipole resonance", which is a resonance peak 
in the cross-section as a function of energy. The dipole state is so 
named because it has a large electric dipole matrix element with the 
ground state. As a result of this last fact, these states lead to strong 
photon absorption cross-sections which are orders of magnitude larger 
than those of other multipoles. 

If nuclear recoil is ignored, the cross-section for electron 
scattering can be written in the general form [Z] : 



dg 
d a) 



n g 2 r 

4 k L 



2kk' cos (e/ 2) f L (q) 



2kk ' . 2 



sin 2 (0/ 2) [ q 2 + 2kk' cos 2 (0/ 2) ] f ^ (q) J 



where q is the four -momentum transfer and 0 is the angle between the 
incoming and final electron momenta, f and f are the longitudinal 

Lj t 

and transverse form factors. Regardless of the model used, the cross- 
section will have this form and will depend on these two form factors. 

If nuclear recoil is taken into account, the cross-section is simply 
multiplied by a factor 

1 

, k 7 - k cos 0 
1 + 

The transverse form factor can be easily separated out, since for 
o 

0 = 180 , the coefficient of the longitudinal form factor goes to zero. 
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Figure 1 shows the experimental results for the square of the 

1 2 

transverse form factor for electron scattering from C as a function 

of momentum transfer, along with the predictions of two collective 

models (the Goldhaber - Teller and Steinwedel - Jensen models) and two 

particle -hole shell models (the Brown and Gillet - Melkanoff models). 

Figure 2 shows the experimental results for the transverse form factor 

for 0 16 and the predictions of the Brown and Goldhaber - Teller models. 

12 1 6 

For both C and 0 , the observed behavior is more closely predicted 

by the shell models than by either of these two collective models, Lewis 
and Walecka assert that this is conclusive proof that the collective mode 
description of these excited states is not adequate, and that the particle- 
hole shell model description is correct. The purpose of this paper is 
to develop a collective model which can predict the observed behavior 
of the form factor. 

In the Goldhaber - Teller model, the neutrons of the excited nucleus 
arc assumed to move against the protons as a unit, thus giving the 
effect of an harmonic oscillator, with the restoring force being due to 
changes in the overlap of the protons and the neutrons [2]. The 
interaction Hamiltonian is assumed to be 

H = j P 2 / u + X - uwQ 2 , 

where w is the frequency of oscillation, u = NZM / A is the reduced 
mass of the system, Q is the coordinate representing the separation 
of the center of mass of the proton system from the center of mass of 
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the neutron system, and P is the canonical momentum associated with 
the coordinate Q. P and Q are quantized and the displacements are 
assumed to be small. The charge and current densities are found from 
the motion of the protons. 

In the shell model theories, a single nucleon is assumed to be 
excited from a closed shell by the interaction with the incoming electron. 
This nucleon then interacts with the resulting n hole M left in the shell it 
previously occupied. The nucleon and the hole are treated as point 
particles. The excited state of the nucleus is assumed to be a linear 
combination of the single-particle excitations, and the charge, current, 
and magnetization matrix elements between the single -particle states 
are found. 

In this paper, a collective mode is assumed and a semiclassical 
treatment is used to find the matrix elements of the Hamiltonian de- 
scribing the interaction between the electron and the nucleus. The 
cross-section is then found to depend on the two expected form factors. 
These form factors are then calculated for a nucleus of constant density 
and for a nucleus whose charge distribution is the Fermi distribution. 
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2. The Model. 



The collective inode assumed is an harmonic oscillation of the 
protons and neutrons in the Q -direction about some axis, so that the 
charge density of the excited part of the nucleus is: 



p ( r ) = Z 7 e f ( r ) ( 1 + a cos© e 1 ^ 1 ) 
where Z ‘ is the number of protons involved in the excitation, and 
Z 7 e f(r) is the ground state charge density of this portion of the nucleus. 
It is assumed that the neutrons' motion is equal and opposite to that of 
the protons, so that the net angular momentum is zero. At this point, 
no specification is made as to whether the entire nucleus or a portion 
only is excited, although it may later be shown that in order to get the 
desired results for the transverse form factor, it is necessary to assume 
that only an outer core is excited, with an inner core of nucleons re- 
maining inert, uu , the frequency of the oscillation, is equal to the 
energy transferred to the nucleus by the interacting electron. The 
system of units employed in this paper is that in which c = “h = 1 . 

From charge conservation, 



V • 



J = 



5p 

at 



The assumed mode has 



J = J e 

e e 



where e is the unit vector in the 0 direction. 

0 



Thus, the charge 



conservation equation is: 



r sin0 d 0 



— — (J sin0) = -iuoaZ^ f(r) cos 0 e 



i a) t 



This leads to: 

J(r) = iuoaZ^ r f(r) sin0 e la)t 

The calculation of the matrix elements H . is done semi-classically. 

fi 

The matrix element is [5] 



H = f J (x) A (x) d 3 x 
fl J |JL |J. 



It is assumed that the nucleus is scattered by the electric field of the 
electron, which can be written as 



a / x 4 Ilie iq * x - 
A (x) = — r e — - u y u 

» q 2 V ‘ ' 



where q is the four -momentum transfer and u„ and u. are the final 
— f l 

and initial plane-wave electron states [5], V is the normalization 

volume for these states. It is further assumed that J (x) for the 

M* 

nucleus can be written as just the classical charge and current densities 
for the model. Thus: 



H 



fi 



,3 4 nie i q * x , - .3 

4 ” - — ~ - J (x) u y u. 

|JL i pi 1 



= Fj (x) A (x) d~ x = ' * * : - w e* — — J (x) u v u. d“ x 

J U LL 4. 



q v 



But J (x) is just j (x) e la>t , and 



iq • x iq * x -iout 

e — — = e e 
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so what must be calculated is the Fourier transform of the current 



density: 



j 



(q) ■ f 



1 q 



x . , > ,3 

J (x) d X . 
M- 
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3. Calculation of Fourier Transforms. 

There will be three normal modes for this model: one for which 

the axis of vibration is parallel to q (the longitudinal mode) and two for 
which the axes of vibration are perpendicular to q (the transverse modes). 
For the longitudinal mode, let the z-axis denote the axis of vibration. 

Then: 



Jx (q> = C 



3 2 iqrcos 0 

r f(r) sin 0 cos0 cos0 e drd0d0 



= 0 



Likewise: 



o 3 2 

j (q) = C r f(r) sin 0 cos 0 sin0 e 

y j 



iqrcos 0 



drd0 d0 



= 0 



Thus, for the longitudinal mode, both transverse terms are zero. 
However: 

j (q) = - C r^ f(r) sin^0 e ^ rcos ® drd 0d0 

z 



8nc 



J f(r) [sinqr - qr cos (qr) ] dr 



= 4IliaZ / ea) f (q) 

L 



where 



f L,( q ) = J ^ sinqr " qr cos ^ qr ^ ^ dr ‘ 



The fourth component of j (q) is: 



. , i q • x 3 P . , v 

p (x) e — — d X = ia j f(r) 



2 . iqrcos 0 

r sin0 cos 0 e drd0d0 



= - 4 naZ'eq f (q) . 
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Thus, for the longitudinal mode, 



j (q) 411 iaZ'eiu f (q) (0, 0, 1, -) , 

g 1^ CL 



For the transverse modes, the calculations are done in the coordi- 
nate system in winch the z-axis is the axis of vibration and the y-axis 
is in the direction of q . For one of these modes: 



, , v 2 lqr sin0 sin0 , 

) p ( q ) ia f(r)r s m 0 cos 0 e drd0d0 = 0 



J x (q) 



^ \ r . K 3 2 iqr sin0 sin0 

C f(r)r sin 0 cos0 cos 0 e drd0d0 



Jy(q) 



c 

0 



r , v 3 2 iqr sin0 sin0 

f(r)r sin 0 cos 0 sin0 e drd0d0 



(q) 



C 



r . . , 3 .3 iqr sin 0 sin0 

f(r)r sin 0 e drd0d0 



who re 



q n c 



| f (r) [ qr - sin (qr) ] dr = - 2 nia a> (q) 



f^(q) = -h J f (r) [ sin (qr) - qr ] dr 



For the modes in which the axis of vibration is perpendicular to q , 
the only nonzero Fourier transforms are those in the directions of 
the axes of vibration. The overall current density Fourier transform 
can be written: 
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j p (q) = 4nu»z'e [- i« pl f t (q). - 5 « p2 f t (q). 6 P3 V q) - 



! 6 P3 f L (q) ] = 4 n.a„z'e f P ( q ) 



where P is the polarization axis for the excited nucleus (1, 2, or 3, 
with the 3-axis being the axis parallel to q ) . 
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4. Calculation of Cross-Section. 



The interaction Hamiltonian is now: 



H 



4 n ie 

q 2 v 



j (q) u f y 



M- 



U. 



.An 2 -7' 2 

1 6 n a Z e at 
2 

q V 



f P (q) u v u 

p, I p, 1 



The cross- section depends on the square of the Hamiltonian, 
averaged over initial electron spin states and summed over final 
electron spins and nuclear polarizations. This is: 




1 

6 



L 

i, f, P 



( 



1 6 n a Z 



2 

e ui 



2 

q v 



)' 



f (q) u v u | 

[1 I }jL 1 



Th< summations lead to a trace calculation, as follows: 



r P ❖ P * ^ 

Tj ( f (q) u y u u . y Y . u r i ) ~ 
rT ^Vu f u i i v 4 f v J 



i, f. P 



Tr \( u. y u. u. y Y. u f ^ P f P 
LV t j, i i v 4 f y v a J 



Let us now examine the quantity in parentheses. If the y and 
Y^ were reversed, this would just be 



Tr ( u y u u - Y u r ) > 

V f Li i i v f y 



which is [ 5 ] : 
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E E 



[kk' + k k ' - (k *k / +m 2 )6 1 

l-« ' U . . \ 1 M ' ' . . \ J 



M,v 



From the commutation relations for the gamma matrices: 



Tr ( u y u. u. y y u ) = 

f 'g, l i ’ v 4 f 



1 



[ k k 7 +k k / 
E.E f m, v 






( k * k + m ) (6 -26.6.) 

— — |j,v p,4 v 4 



The form factor matrix is: 



f = z f2 (q) 

4 t M 



1 ,2 



10 0 0 
0 10 0 
0 0 0 0 
0 0 0 0 



+ f L (q) 



* 4 f , «*> \„ + f L «» A vu, 



0 0 0 



0 0 0 



0 0 1 



iq/uj 



2 2 

0 0 -iq/u) q / a» 



There are two traces which must be evaluated. These are; 



The transverse coefficient: 

Tr (B A ) = 4 sin 2 (9/2) f 1 

Vn n,v' L 



2kk / 2 , .„ 1 

— cos (9/2) J 



The longitudinal coefficient: 



Tr (B A ) 

Vp, nv 



2q 



cos (9/2) 



where 0 is the angle between the incoming electron momentum k and 
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the outgoing electron momentum k'. The square of the Hamiltonian, 
averaged and summed over states, is: 



ltr , 2 r 1 6 n " az'e N 2 2 1 

|H | = ^ -5 ) <u 



q^V 



) " TTEP { 2kk ' sin (e/2 > 



r , 2kk / 2 , /ov 1 2 , x 

X [ 1 + — cos (0/2) J f t (q) 



+ 4kk ' q^ 

2 2 
q u) 



COS Z (e/2) (q) | . 



The amplitude a of the oscillation can be calculated (classically) 
from energy considerations. At the time during the oscillation such 
that the charge density of the nucleus is its ground-state configuration 
( a) t (n + ~ ) n ) , all the energy transferred to the nucleus is in the 
kinetic energy of the neutrons and protons. The velocity distribution 
of the excited protons is 

V (r) = "j (r) / Z'ef(r) . 

This ratio is v(r) = ^ uoar sin0 . The mass density of the excited 
protons is Z / M f(r). Thus, the kinetic energy distribution of the 
excited protons is: 

MZ* f(r) v^(r) = ~ Z / M f(r) uu a r sin 0 . 

c o 

Since half of the transferred energy will be in the protons (with the 
other half being in the kinetic energy of the neutrons): 
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(jo MZ'uo a f r/ x 2 .2 2 . 

= - f(r) r sm Q r sin0 d 0 dr d0 

2 8 J 



This gives 



2n Z'Muol 



where I = 



= j f(r » 



r 4 dr. 



The cross-section is, from the well-known Golden Rule of Fermi, 



d a = L 

final 

states 



2n 



I Hi 2 8 (E - E. ) — — 
1 fi 1 u v f l ' flux 



The flux is just k/ VE. . The summation over final states is just 

3 / c 

an integral over d k for the case at hand. The O-function is just 



6 (k - k' - u) ) . Thus: 



da 



VE 

V i 



( 2 n) 



2 k 



| H | 2 6 (k - k' - u>) k 7 2 dk # dQ 



v 2 e. k /2 

I I 



(2n) k 



k - k = ua 



Finally, putting |H| into this equation, we obtain: 



. 2 , 

da _ Of k 

d w 4, 

q k 



— ^ I 2kk ' sin 2 (0/2) f (" 1 + 2k ^ cos 2 (0/Z) 

n 2 Z'MI 1 L q 



4kk ' q 4 

_ 2 2 ~ 

q (jo 



2 2 
f cos 



(e/2) j 
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Note that these form factors must be multiplied by a constant 



term 



zn z z'm i 



in order to compare them with the form factors shown in Figures 1 
and 2. 
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5. Uniform Shell Model. 



If the nuclear charge density is assumed to be uniform, with the 
nucleons within a shell of inner radius r^ and outer radius r^ being 
excited by the interaction, we have: 

f L ( q ) = -ij ["sin (qr^ ) - qr^ cosfqr^) - [sinfqi^) - qr ^ cosfqr^jj 

q 

f t (q) = -h r s in(qr 2 ) - qr 2 - [sin(qr 1 ) - qi^] J . 

q 

Figure 3 shows the results for the square of the transverse form 

factor versus momentum transfer for 0^ , with the assumption that the 

entire nucleus is excited (r^ = 0) , and the same form factor with the 

4 

assumption that an inner core consisting of a He nucleus remains 
inert. The figures used for the outer radii for these calculations were 
3.41 fermi for 0^ and 2.08 fermi for He^ [ 3 ] . 

If one assumes an even thinner shell is excited, the form factor is 
reduced by orders of magnitude and falls off to zero with increasing q 
even faster. Figure 4 shows the results for a shell of thickness 0. 1 
fermi. No matter how thin the shell, the uniform shell model cannot 
give the observed form factor behavior with q for a reasonable outer 
radius. 



22 



6. The Fermi Distribution Model. 



Now let us assume that the charge density function is just the 
Fermi distribution: 



f(r] 



1 



1 



I ,." 1 " 11 lte a(r - b) 



The integral can be done by breaking it up into two parts and writing 
each part as an infinite series, as follows: 



f Cl) 



f (r) [ sin (qr) - qr ] dr 



r r a 


'• CO 


+ 




L J n %} 


A “* 



f(r) 



X [ sin (qr ) - qr ] dr 



For r < a : 



1 



+ e 



a(r - a) 



£ (-1)" e "“< r -a) 



For r > a . 



1 



1 l e 



(r -a) 



» ,. 0 nH e -na-(r-a) 



Tin. r insult for the form factor is: 



C (q) = 






2 2 

n 



6cv 



2 J 



2 sin 






- n a a 



2 cos q a 



- L ( - l ) 
l 



n j" n 



2 2 

O' 



2 2 2 
q + n a 
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For excitation of a shell, the form factor is 



* t (q) = f. (q) - f t (q) 



where a is the outer radius of the shell and b is the inner radius. 

1 6 

The square of the transverse form factor for 0 is shown in Figure 5, 

using a = 3. 64 and these combinations for b and a: b = 3. 0 , 

a = 1; b = 3. 6 , a = 1 ; and b = 3. 0 , a = 20 . Again, the thinner 

the shell, the smaller the form factor and the faster it drops off with 

increasing q. For large a, the square of the form factor is smaller 

and does not drop off as fast. In fact, for very large a, one gets just 

2 

the uniform shell. Note that the f curve for a = 20 actually crosses 
the a - 1 curve. This leads one to suspect that if one assumes an outer 
boundary for the excited shell with, say, a = 1 and an inner boundary 
with very large a, for a thin enough shell one could get a curve for 

i i 2 

| f^ I which has the desired minimum at about q = 0. 3 . 

Figure 6 shows the results for three trials of this composite model. 
For a thick shell (a = 3. 64 and b = 2. 0) with a = 1 , the form factor is 
larger than for the same shell with a = 10. The form factor becomes 
larger for a thinner shell. However, for the combinations of a, b, and a 
used here, the form factor does not get the desired behavior. 
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7. Conclusions. 



None of the assumed charge distributions gives the desired behavior 
for the form factor. This model gives somewhat the same results as 
the Goldhaber - Teller model; that is, the function | f (q) | is a 
motonically decreasing function of q within the range of interest. 

Further work is necessary on the composite model before this model 
can be completely discarded, however. For instance, if one assumes 
the nucleus as a whole has a Fermi charge distribution and that the inert 
core is a uniform shell with a radius equal to the radius of the half- 
height of the Fermi distribution, what one has is essentially a classical 
treatment of the part icle -hole theory. Figure 7 illustrates the 
"particle" and the "hole" associated with this model. This idea should 
be pursued further. There is a strong possibility that it will give re- 
sults similar to those of the successful Brown particle-hole theory. 
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I 




Figure 1. Square of the transverse fom factor for 
the giant clipole resonance for Carbon 12. 
{ Lewis and Walecka, 19 6 )',) 
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I'lCure 2. Square o±‘ the transver.ee rom factor ror 
the giant dipole resonance i'or Cec'Ccn lo. 
( Lewis, ly63; 
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u u.2 u.U 0.6 u.d l.U 1.2 l.u 



Momentum Transfer q j 

figure 3. Square of trie transverse fron factor for 
Oxy-on 16, entire nucleus excited, and 
outer shell only excited, uniform shell 
model. 
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1*1 Cure U, Square of the transverse forr. factor for 
Oxygen lo, very thin ^.1 f; shell of 
constant density*. 
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Momentum transfer q 



ijl^arc q. Square of the transverse form factor for 
Oxygen 16, Feimi distribution charge 
density, o<— 1 and oi — '^O • Tne dashed 
line represents loo times the square of 
the transverse form factor for a very 
thin shell ^*olj. fermi ,1 with o(= 1 • 
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Figure 6, Square of the transverse forr factor for 
Qr® for tte composite model. 
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Figure Y. Classical particle- nole analogue. 
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lective models. The purpose of this paper is to develop a collective 
model which could possibly give this dip in the form factor. Further 
work is necessary to determine if this model will actually give this result. 
A semiclassical treatment is used to find the matrix elements for the 
assumed model. The transverse form factor is then calculated for two 
nuclear charge density functions -- constant density and the Fermi 
distribution. 
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